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Abstract. This paper describes the partial wave expansion and inte- 
gral representation of Bessel beams in free space and in the presence 
of dispersion. The expansion of the Bessel beam wavepacket with con- 
stant spectrum is obtained as well. Furthermore, the sum of a triple 
Legendre polynomial product of same order but different argument fol- 
lows naturally from the partial wave expansion. The integration of all 
Bessel beams over all conical angles is shown to have a simple series 
representation, which confirms the equivalence between the results for 
both expansion and integral representation. 



1. Introduction 

The Bessel beam is a mathematical construct that is a solution to wave 
equation. Bessel beams have attracted attention ever since Durnin [1] pub- 
lished his paper in 1987. The attention is due to the beam's intensity profile 
that has a sharp peak along the axis of propagation. The electromagnetic 
Bessel beam also travels with a phase velocity that is higher than speed 
of light. Numerous authors have used a superposition of Bessel beams in 
order to form a wavepacket (for example see [21 13 H] ) • Their choice of using 
an exponential spectrum has allowed them to create a wavepacket that has 
been called an X-wave. The X-wave however is not a physically possible 
solution because it carries infinite energy. There is also interest in the su- 
perluminal property of X-waves. Two papers have been recently published 
that clearly show how two schools of thought concerning the X-waves have 
opposing views. One side states that X-waves are physically realizable, even 
superluminally [5], whereas the other side states that superluminal effects of 
X-waves are simply phase related and that superluminal information trans- 
fer via X-waves is not possible [6]. The author of this paper shares the latter 
view. Nevertheless, the mathematics behind Bessel beams, which are build- 
ing blocks for wavepackets, is very rich as shown in this paper. The mathe- 
matical tools developed here may help create a proper wavepacket that has 
finite energy and is therefore physically realizable. Such a wavepacket may 
still be of interest because Bessel beam itself has reconstructive properties 
not recognized in other types of beams [7j. In this paper I will show and con- 
firm two representations of Bessel beam, namely a partial wave expansion 
and an integral representation. The partial wave expansion of Bessel beam 
has been reported before, as will be explained later, but the proof presented 
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here is new and based on a physical argument. The integral representation 
however is to my knowledge a new result. Furthermore I will show that a 
simple superposition of Bessel beams with a constant spectrum is supris- 
ingly related to the infinite series of triple Legendre polynomial product. In 
the process I will show a Bessel beam identity that is simple and also new. 
I will also briefly describe the above two representations when dispersion is 
included into wave equation. 

2. Mathematical Introduction 

A zeroth order Bessel beam that propagates in the z direction can be 
described by 

(2.1) ^(/o, z, t) = exjp{ikzZ — iujt)jQ{kpp) 

where p is the transversal distance from the propagation axis z. 
Solution (12. ip satisfies the wave equation 

,2,2, V^.4f 

where 

(2.3) kz = - cos(e) 

c 

(2.4) kp = -smie), 

c 

which is a wave vector construct used by Salo et al. [3]. 

The angle 6 is the axicon angle of the wave. The wave vector k = 
{kz,kp,k^) is restricted to lie on a cone described by vertex angle 6. The 
Bessel beam itself constitutes the sum of all wavevectors on this cone (i.e. 
sum over all A:^). For simplicity I shall let c = 1. 

I will first show that the partial wave expansion of a Bessel beam is 

exp[iu) cos{9)z — tujtjJoluj sin{6) p] 

(2.5) = V 2z"(n + -)Pn[cos{9)]Pn[ J-— \Uu:^/z^ + p^) exp[-^t^t], 

^0 2 ^Jz^ + P^ 

where Pn is Legendre polynomial of order n and jn is a spherical Bessel 
function of same order. 

I will then show that the integral representation of Bessel beam is 

exp[2a; cos{6)z — vjjt]jQ[oj s\n{9)p] 

1 f°° izX 

(2.6) =-/ jo(^)exp[ ] exp(-^a;t)dA, 

71" J-oo V2:2 _^ p2 

where jo is spherical Bessel function of order zero. R is given by 



(2.7) R = V A2 + w2(z2 + p-i) - 2Alj V^2 + p2 cos(0). 
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The Bessel beam wavepacket is a superposition of Bessel beams: 

/oo 
S{u:) exp[^a; cos(6')z — xwt] Jo[a; sm{6) p]duj . 
-oo 

If S{uj) = 1,1 will show that 

exp[iuj cos{6)z — iujt]jQ[uJsiyi{9)p]duj 



oo 

(2.9) = ^==^27r(n + -)P„[cos(e)]P„[cos(7?)]P„[cos(7)]. 

where cos(??)= . ^ and cos (7) 



The integral in (j2.9p can be evaluated analytically, and from this I will 
show the result for the following triple Legendre polynomial sum 



(2.10) 



^(2n + l)P„[cos(e)]P„[cos(7?)]P„[cos(7)] 

n=0 

_ 2 

\/ sin^(7?) sin^(0) — [cos(r/) cos(0) — cos(7)]^ 



if sin(r/) sin(6') >| cos(r/) cos(6') — cos(7) |, otherwise the sum equals zero. 

The result in (I2.10p is elegant, short and simply obtained. Similar result 
has been found by Baranov [TO] by using different, more complicated and 
lenghty method. My result coincides with that of Baranov jl9] and I have 
numerically verified it as well. 

I will finally evaluate the following integral: 



(2.11) / Jo[ujrsm^{e)]exp[iujrcos'^{9)]d[cos{e)] = S22z''jn[ujr) 

■^-1 n=0 



For this I will use both the partial wave expansion and the integral rep- 
resentation of Bessel beams. In both cases I get the same result and this 
verifies the equivalence between (|2.5p and (|2.6p . The result itself constitutes 
a Bessel beam identity which to my knowledge is new. 

Marston [8] gives and numerically verifies the Bessel beam partial wave 
expansion in (12. 5p . Mitri [9] shows the Bessel beam partial wave expansion 
directly from the Gegenbauer result published in 1899. I can however find 
the exact form of the expansion as in (12. Sp in the work of Stratton [10^ 
page 413]. The proof of (12. Sp in this paper is new, and is described from a 
physical perspective. The integral representation in (12. 6p is to my knowledge 
new. 
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3. Bessel Beam Partial Wave Expansion 
Let me write a superposition of Bessel beams with varying conical angle 

6: 

(3.1) j -B[cos(0)] exp[zLi; cos(0)z — zwt] Jo[t^ sin(6')/9](i[cos(0)], 

where B[cos{9)] is angular spectrum. The angular spectrum can then be 
expanded in terms of Legendre polynomials, namely 

oo 

(3.2) B[cos{e)] = J2 anPn[cos{e)], 

n=0 

where Pn is Legendre polynomial of order n and coefficients o„ are constant. 

If I choose the spectrum to be a delta function that eliminates all Bessel 
beams from the superposition except the one with conical angle 

(3.3) B[cos{e)] = 6[cos{e) - cos(eo)], 

then one has to let a„ = (n + ^)P„[cos(0o)]) because the standard identity 
for Legendre polynomials is: 

oo ^ 

(3.4) ^(n + -)Pn[cos{eo)]Pn[cos{9)] = 6[cos{e) - cos{9o)]. 

n=0 

I let a = cos(0) and oq = cos(0o)- At t = and with delta function 
angular spectrum, (13. Ih becomes: 

(3.5) / 'S~]{n + -)Pn{ao)Pn{a)ex.p{iuaz)Jo{u\/l^^a^p)da. 

•^-^ n=0 

The expression in ([33]) is equal to the expression for Bessel beam at t = 0: 



exp{iu}aoz)Joiuj\J 1 - a^p) 
(3.6) = / y^(n + ■^)Pn{ao)Pn{a) exp{iu}az)Jo{u}y/l - a'^p)da. 

-'-^ n=0 

The right hand side of p.6p is very usefull because the integral in (j3.6p 
can be evaluated, (see [ini [HI HI [l3l [Hj ) . I should mention that Neves et 
al. [11] state that they do not know of any other report of this evaluation 
result. However Cregg and Svendlindh [12], Koumandos [13] and Dodonov 
|14j clearly show that the result follows from publications of Gegenbauer at 
the end of 19th century. Nevertheless, Cregg and Svendlindh [12J state that 
the result is not generally well known. Stratton [iHl page 411] evaluates the 
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integral in (13. 7|) in the exact form as well as authors in [HI [121 [HI E! ■ 

Pn{a) exjp{iujaz)jQ{io\/l — oPp)da 

1 

(3.7) = 2^"P„[^^=]i„(.;7^^). 



Equation (|3.6p therefore becomes the partial wave expansion of Bessel beam 
written below for any time t: 



exp(za;ao2 — iojt)jQ{ujy 1 — a^p) 

oo ^ 

(3.8) = V 2i"(n + -)P„(ao)^n[^=^b„(u;V^2 + ^2) exp(-zwi). 

n=0 

4. Bessel Beam Integral Representation 

I can simplify the Bessel beam partial wave expansion in ()2.5p at i = if 
I let j3= I I and p = + obtain 

/ °° 1 

(4.1) exp(^L^ao2) Jo(wV^ " "o-") = ^'"^'^ + -)P„(ao)Pn(/3)in(^). 

?i=0 

The Fourier transform relationship between Legendre polynomial -Pn(/5) 
and the spherical Bessel function jniX) is given by Mehrem [15j : 

(4.2) p„(/3) = ^^/ j„(A)exp(^/3A)dA. 

J-oo 

I can therefore further rewrite equation (14. ip as 
exp(2(jjao-2) Jol'^Y^l — (^oP) 

(4.3) =-/ Y.{n + -)Pn{ao)jn{X)jn{p)exp{i/3X)dX. 

The spherical Bessel function is defined as jn(A) = ^/^Jn+lW^ where 
J^_,_i (A) is the Bessel function of half integer order. With this, (j4.3p becomes 

exp{^u;aoz)Jo{uJ^J 1 - a^p) 

(4-4) = y ^|E(^+2)^"(«°)^% l^Uxp(z/3A)dA. 

The equation (|4.4p is very useful because the sum within is evaluated by 
Hochstadt [HI page 221]), 

(4.5) 5](n + ^)P„(ao)^^i^= = -Jo{R), 
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where R is given by Hochstadt [Ml page 219] 
R = + /x2 - 2\^iao 



(4.6) = y A2 + ^^2(^2 + p2) _ 2XoJ^/^Tf^cos{eo). 

I have therefore reduced equation (j4.4|) to the integral representation of 
Bessel beams written below for any time t. 

exp{iujaoz — iu:t)jQ{ijj\J 1 — a^p) 

1 /■°° izX 

(4.7) =-/ jo(^)exp[^==]exp(-2a;t)dA. 

71" J-oo Vz2 + p2 

5. Expansion of Bessel Beam Wavepacket with Constant 

Spectrum 

Bessel beam wavepacket with constant spectrum can be expanded into an 
infinite series of product of three Legendre polynomials. For this purpose let 
us integrate the partial wave expansion from (|2.5p over variable u in whole 
frequency domain: 

(5.1) 

exp[zcL! cos{6)z — iojt]jQ[oj s\n{9)p]dijj 



oo 
oo 



^ 1 />oo 

V 2z"(n + )p„[cos(0)]P„[^^] / U^yjz^ + p^) eM-^^t)du. 



n=0 

I also make the following variable substitutions 
(5.2) cos(r?)= 



(5.3) u 



Vi2^ 



The right hand side of (j5.ip then becomes 
(5.4) 



^== ^ 2z"(n + -)P„[cos(0)]P„[cos(7?)] / M exp[-^L=]d/i. 



I recognize the integral inside (j5.4p as Fourier transform of Legendre poly- 
nomial namely 

(5-5) Pn\ I I J = , \„ / in(/^)exp[— 
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I let cos (7) = , * and therefore (|5.1|) becomes the desired expansion 



exp[za; cos(6')z — iujt]jQ[uj sin{6)p\duj 

00 

(5.6) = ^==^2^(n + -)P„[cos(e)]P„[cos(7?)]P,,[cos(7)]. 

V + „=o ^ 

The left hand side of (15. 6p is a superposition of all Bessel beams with 
frequencies uj from —00 to 00, and whose frequency spectrum is /(w) = 1. 
The integral on the left hand side of (I5.6p can be evaluated analytically as 
well and the result is: 



(5.7) 



exp[iu} cos{9)z — iujt]jQ[ijj sm{9) p]doj 

00 

2 

Vsin^(%2 _ _cos(0)z]2' 



The integral result in equation (j5.7p is valid if | t — cos{9)z \< sm{6)p, 
otherwise the integral is equal to zero. This means that the expansion in 
(|5.6p vanishes if sin(7;) sin(0) <| cos (r/) cos (0) — cos(7) |. I therefore have 
evaluated the sum of the triple Legendre polynomial product as well: 

00 ^ 

^ 27r(n + -)P„[cos(e)]P„[cos(r?)]P„[cos(7)] 

n=0 



(5.8) - 



y^sm^{e)p^ - [t-cos{9)z]^' 
If I use the definitions for cos(7) and cos(r7), I finally get 



00 ^ 

J2 Mn + -)Pn[cos(^)]P„[cos(r?)]P„[cos(7)] 

(5.9) 



n=0 

2 



Y^sin2(r/) sm'^{9) — [cos(??) cos(6') — cos(7)]2 

The result of (j5.9p is valid if sin(r7) sin(^) >| cos(??) cos(6')— cos(7) |, otherwise 
the sum is equal to zero. The expansion of Bessel beam wavepacket with 
constant spectrum is a new result as well. 



6. Bessel Beam Identity 

In order to show that the two representations in (j2.5p and (j2.6p of Bessel 
beam are equivalent, I shall consider the following case: I will evaluate both 
representations along points of the cone only at t = and I will integrate 
the results over all possible cones. Let me start with expansion first. At a 
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point on the cone I have cos(r/) = cos{9). I let r = \J ^ and I write: 
J{i\ujir sin^(0)] exp[zwr cos^(0)](i[cos(0)] 

n=0 
-1 oo 

3.1) = / J^2z"(n + -)j„(a;r)P„[cos(^)]P„[cos(0)]d[cos(^)]. 

•^-1 n=0 

If I exploit the orthogonality of Legendre polynomials I get 



oo „i 

^2i^{n+-)Uu;r) P„[cos(0)]P„ 



cos(0)]ci[cos(0)] 



n=0 

oo 



^2z"(n + i)j„(a;r)-^ 
^0 2 + 2) 



00 

(6.2) = J^2z"i„(^^r). 

n=0 

Similarly, let me proceed wih integral representation: 

'1 



Jo[c<jrsin (^)] exp[zLi;r cos {9)]d[cos 



-1 



/I 1 /"OO x 
- / io(i?)exp[^==](iA(i[cos(0)] 



-1 " J-oo " Y^Z^ + p2 

- / io(^)exp[?Acos(7)]dAd[cos(6l)]. 

-1 ^ J-OO 

The well known plane wave partial wave expansion namely [17j 

00 ^ 

(6.4) exp[zAcos(7)] = J^(n + -)z"j„(A)P„[cos(7)] 

n=0 

is used to rewrite the last line of (16.31) as: 



/' ^ io(i?) |f;(n + ^)z"i„(A)P„[cos(7)]| (iA(i[cos(0) 

1 1 />oo r /•oo 

(6.5) =-Y^{n + -Y Jn(A) / io(^)^'n[cos(7)](i[cos(e)]|>dA. 



The inner integral of the right hand side of equation (|6.5p can be evaluated 
THl and therefore (16.51) reduces to: 



(6.6) -Y.(n + -y J ^j„(A)|^J„+i(A)J„+i(^r)|dA. 
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Further simplification leads to 

Y.^ri + -)ejn{ujr) \ d\ 

n=0 ^ ^ 

-1 n l-OO 

(6.7) = 5^(n + -)z"j„(u;r)- / UXfdX. 

n=0 ^ ^ ''-^ 

The last integral in (j6.7p can be evaluated as well: 

/■oo 

(6.8) / j„(A)2dA 



2n+l 

My final result therefore is 

„1 oo 

(6.9) / Jo[<:i^rsin^(0)] exp[za;r cos^ (0)]d[cos(0)] = J^2^"j„^a;r) 

•^-^ n=0 



which is the same result as in (|6.2|) . The equivalence between the partial 
wave expansion and integral representation of Bessel beam is therefore ver- 
ified. The result in (|6.9p constitutes a new Bessel beam identity. 



7. Partial Wave Expansion and Integral Representation with 

Dispersion 

The importance of the partial wave expansion of Bessel beam becomes 
more obvious when dispersion is allowed into the wave equation. Dispersion 
arises when the index of refraction depends upon frequency. In this case the 
Bessel beam partial wave expansion becomes 

(7.1) 

e-xp[in{(jj)uj cos{6) z — iojt\jQ[n{u))uj sm{0)p\ 

oo ^ 

= V 2z"(n + -)P„[cos(g)]Pn[ / ]jn[n{uj)oj^z^ + p^] exp[-zwt], 

where n{uj) is frequency dependent index of refraction. 

The index of refraction in the left hand side of (|7.ip appears inside the 
Bessel function as well as inside the exponential, whereas n{(jj) appears only 
inside the spherical Bessel function in the right hand side of (j7.ip . If a 
certain form of n{ijj) allows the summation in (j7.ip to be terminated at a 
certain integer value, then this may allow for better approximations. In the 
case of integral representation n(a;) is present inside the variable R only: 



(7.2) R = ^A2 + n2(6<j)a;2(z2 + ^2) _ 2An(w)a;v^ + p2 cos(0). 
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8. Conclusion 

In this paper, I gave a new proof for partial wave expansion of a Bessel 
beam, from a physical perspective. The partial wave expansion has been 
shown by Stratton [10^ page 413] in a more general form. To my knowledge 
the integral representation of Bessel beam however is not found in literature. 
Furthermore the equivalency between the Bessel beam wavepacket with con- 
stant spectrum and the triple Legendre polynomial series as shown in section 
5 is a new result as well. The similar sum has been found by Baranov [19] by 
using different more complicated method and our result coincides. The final 
result of section 6 further shows the equivalence between the partial wave 
expansion and integral representation of Bessel beams, which further con- 
firms the integral representation. I have further developed the mathematical 
structure of Bessel beams with the hope that the experimental results in 
Bessel beam optics can be described more accurately and the controversy of 
superluminal propagation resolved. 
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